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ABSTRACT 

Anomalous  cloud  lines  are  frequently  seen  in  satellite  images 
as  curvilinear  features.  These  cloud  lines  or  "ship  tracks"  are 
likely  due  to  products  of  ship-produced  hot  exhaust  gases  that  are 
expelled  into  the  atmosphere,  increasing  the  aerosol  concentration 
in  the  ship  track  plume.  NOAA  9  and  10  AVHRR  data  are  sensitive  to 
cloud  droplet  size  and  show  the  ship  tracks  as  increases  in 
radiance  due  to  reflectance.  Twenty  eight  NOAA  9/10  satellite 
passes  are  analyzed.  Twenty  two  of  the  passes  are  found  to  contain 
a  total  of  316  ship  tracks  which  is  significantly  more  than  that 
expected  by  earlier  ship  track  studies.  An  existing  ship  track 
detection  algorithm  is  used  to  conduct  a  statistical  comparison  of 
ship  track  and  non-ship  track,  or  ambient  pixel  reflectance  of  the 
NOAA  9  and  10  AVHRR  channels  1  (0.63  Mm),  3  (3.7  /xm)  ,  and 
4  (11  /xm) .  The  results  of  the  statistical  analysis  confirm,  as 
found  in  previous  studies,  that  the  ship  track  pixels  displayed  a 
significant  increase  in  values  for  channels  1  and  3  and  a  very 
slight  increase  for  channel  4. 
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I.   INTRODUCTION 

A.   BACKGROUND 

Much  attention  has  been  given  to  anomalous  cloud  lines 
seen  in  visible  and  near-infrared  satellite  images.  The 
source  of  these  anomalous  cloud  lines  were  first  identified  by 
Conover  (1969)  as  exhaust  from  ships  transiting  at  sea  and  are 
often  referred  to  as  ship  tracks.  He  showed  that  the  most 
likely  cause  of  ship  track  clouds  is  the  increase  in  the 
aerosol  content  produced  by  ship  exhaust.  Ship  exhaust  is  a 
source  of  aerosols  that  act  as  cloud  condensation  nuclei  (CCN) 
that  have  the  effect  of  increasing  the  number  of  cloud 
droplets  while  reducing  the  droplet  size  (Coakley  et  al., 
1987)  . 

Advanced  Very  High  Resolution  Radiometer  (AVHRR)  data  was 
used  in  the  Coakley  et  al.  (1987)  study  with  Channel  1 
at  0.63  /nm,  Channel  3  at  3.7  /xm,  and  Channel  4  at  11  /xm. 
Coakley  et  al.  suggested  that  the  effects  of  aerosols  on  cloud 
reflectivity  may  have  a  much  greater  influence  on  the  earth's 
albedo  than  that  due  to  the  direct  scattering  and  absorption 
by  aerosols  alone.  They  developed  an  algorithm  to  detect  ship 
tracks  automatically  by  comparing  reflectivities  of 
contaminated  clouds  to  those  of  similar  non-contaminated,  or 
ambient  clouds.   The  Coakley  et  al.  algorithm  worked  well  in 


areas  of  uniform  cloud  cover,  but  did  not  perform  well  in 
areas  without  ship  tracks,  cloud  free  areas  or  regions  of 
transition  from  one  cloud  regime  to  another. 

Norehead  (1988)  evaluated  the  Coakley  et  al.  algorithm 
incorporating  improvements  in  its  ability  to  detect  ship 
tracks.  The  modified  algorithm  was  used  in  this  study  to 
statistically  analyze  and  compare  the  radiative  properties  of 
the  ship  tracks  found  in  data  gathered  from  1  to  19  July  1987 
in  the  North  Pacific  Ocean  Basin.  Three  hundred  and  sixteen 
ship  tracks  were  observed  in  this  study.  The  Morehead 
modified  algorithm  failed  to  identify  all  the  ship  tracks  that 
were  found  by  visual  inspection.  It  also  had  a  further  draw 
back  of  identifying  non-ship  track  segments  as  ship  tracks. 
A  manual  filtering  technique  was  developed  for  this  thesis  to 
ensure  that  the  statistical  analysis  was  conducted  only  on 
confirmed  ship  track  segments.  This  process  will  be  described 
in  detail  in  Chapter  II.  An  additional  focus  of  this  study 
was  to  identify  and  plot  the  geographical  location  of  the 
observed  ship  tracks.  This  process  will  be  explained  in 
Chapter  III. 

B.   AEROSOLS 

Aerosols  alter  the  atmosphere's  ability  to  scatter  and 
absorb  solar  radiation.  They  also  affect  the  radiation  budget 
through  their  effect  on  clouds.  The  ability  of  aerosols  to 
effect  clouds  comes  from  their  function  as  CCN  which  serve  as 
anchors  on  which  cloud  droplets  can  form.    Atmospheric 
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particle  size  distributions  cover  a  large  range.  CCN  have  a 
radius  of  about  0.1  urn,  cloud  droplets  have  an  average  radius 
of  10  /xm,  and  large  cloud  droplets  have  a  radius  of  100' s 
of  /Ltm.  While  CCN  are  too  small  on  their  own  to  have  a 
significant  interaction  with  incoming  solar  radiation,  once 
cloud  droplets  form  around  them  their  interaction  with 
incoming  solar  radiation  is  very  much  greater. 

Ships  are  a  source  of  aerosols  in  the  atmosphere.  Exhaust 
expelled  onto  the  atmosphere  is  in  the  form  of  water  vapor, 
gaseous  constituents  and  aerosols.  Aerosols  are  of  the  size 
to  act  as  CCN,  as  discussed  above,  and  boost  the  CCN 
concentration  far  above  the  normal  level.  Twomey  (1968,  1984) 
and  Hindman  et  al.  (1977)  have  described  how  the  concentration 
and  size  distribution  of  CCN  have  a  direct  impact  on  the 
formation  of  cloud  droplets.  The  new  source  of  CCN  from  ship 
exhaust  will  cause  a  change  in  reflective  properties  of 
clouds.  Twomey  and  Cocks  (1982)  show  that  an  increase  in 
droplet  concentration  causes  a  decrease  in  the  size  of  the 
droplet.  The  shift  in  droplet  size,  along  with  an  increase  in 
concentration,  will  result  in  a  significant  increase  in 
reflectance  due  to  scattering. 

C.   CLOUD  MICROPHYSICAL  PROPERTIES  AND  RADIATION 

Electromagnetic  energy  reacts  with  the  atmosphere  either 
by  scattering  or  absorption.  Energy  that  is  absorbed  is 
converted  to  an  increase  in  the  temperature  of  the  absorbing 
body.   The  atmosphere,  having  a  temperature,  can  also  emit 
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electromagnetic  energy.  These  processes  represent  the 
interactions  that  take  place  in  the  atmosphere.  The  source  of 
radiation  reaching  a  satellite  is  then  due  to  either 
reflection  or  emission.  Reflection  results  from  energy 
entering  the  atmosphere  and  undergoing  single  or  multiple 
scattering  events  until  it  is  directed  back  out  the  top  of  the 
atmosphere,  primarily  affecting  short  wave  solar  radiation. 
Energy  which  is  absorbed  is  converted  into  an  increase  in 
temperature  of  the  absorbing  body  which  can  subsequently 
undergo  emission  at  a  multitude  of  wavelengths  including  long 
wave  frequencies. 

1.  Absorption  and  Emission 

The  atmosphere  has  a  direct  impact  on  satellite 
measured  radiance  due  to  absorption  and  emission.  Clouds  are 
composed  primarily  of  liquid  water,  so  absorptive  properties 
of  water  are  a  primary  concern.  Liquid  water  absorbs 
electromagnetic  energy  at  all  wavelengths.  The  absorption  is 
much  greater  at  the  infrared  wavelengths  than  in  the  visible. 
This  leads  to  greater  emission  at  the  longer  infrared 
wavelengths. 

2 .  Scattering 

The  type  of  scattering  interactions  depend  on  the 
cross-sectional  area  of  the  particle  with  which  the  energy  is 
interacting.  When  the  scattering  particle  size  and  the 
wavelength  of  the  incident  energy  are  of  the  same  order,  the 
interaction  is  best  described  by  the  Mie  scattering  theory. 


For  the  particles  in  this  category,  the  scattering  properties 
become  less  dependent  on  the  wavelength  and  more  dependent  on 
particle  size.  Mie  scattering  has  a  greater  percentage  of 
scatter  in  the  forward  direction  than  to  the  side  or  backward 
direction.  As  the  droplet  radius  increases,  forward 
scattering  increases  and  backward  scattering  decreases. 
3.   Reflectance 

In  this  study  we  are  considering  reflectance,  the 
measure  of  energy  reflected,  from  cloud  surfaces.  At  the 
wavelengths  considered,  the  reflectance  depends  most  heavily 
on  Mie  scattering  processes.  A  useful  measure  of  the  degree 
of  scattering  can  be  found  by  examining  the  optical  depth  of 
the  cloud  layer.  Twomey  (1977)  concluded  that  pollution  (or 
ship  stack  exhaust  in  the  case  of  this  study) ,  by  increasing 
cloud  nucleus  concentration,  hence  increasing  the  numbers  of 
cloud  droplets;  leads  to  increased  cloud  optical  depth  which 
increases  cloud  reflectivity.  Hunt  (1972)  described  the  phase 
scattering  characteristics  as  well  as  the  particle  size 
effects  on  reflectance  of  cloud  layers  in  the  visible  and 
infrared  windows.  He  showed  that  the  changes  in  cloud 
radiative  properties  at  the  visible  and  infrared  wavelengths 
are  consistent  with  changes  in  the  size  distribution  of 
particles  such  that,  as  the  radius  of  the  cloud  droplets 
decreases  the  reflectance  increases. 


D.   SHIP  EXHAUST  PRODUCED  CLOUDS 

Ship  tracks  appear  as  curvilinear  features  at  visible  and 
near-infrared  wavelengths  in  satellite  imagery.  They  have  a 
plume  like  nature,  narrow  at  the  source  and  spreading 
horizontally  with  distance  from  the  source.  Ship  tracks  can 
be  several  hundred  kilometers  in  length  and  can  last  for  days. 

It  appears  that  there  are  two  different  types  of  ship 
tracks.  The  first  are  seen  in  the  visible  satellite  imagery 
where  anomalous  cloud  lines  form  in  areas  where  cloud 
formation  is  suppressed.  These  may  be  produced  as  the  ship 
enters  an  area  where  the  atmosphere  is  experiencing  a  CCN 
deficit.  As  the  CCN  from  the  exhaust  gases  mix  with  the  air, 
the  deficit  is  reduced  enough  for  clouds  to  form.  Another 
theory  is  that  the  cloud  free  region  is  in  an  area  where  the 
cloud  droplets  are  of  sufficient  size  to  precipitate  out 
therefore  keeping  the  area  free  of  clouds.  As  the  ship  passes 
through  the  area,  the  addition  of  the  exhaust  gas  CCN  reduces 
the  size  of  the  droplets  so  that  they  are  no  longer  large 
enough  to  precipitate  out  and  thus  the  cloud  forms 
(Albrecht,  1990). 

The  second  type  of  ship  track  can  be  seen  in  the  near- 
infrared  (NIR)  satellite  imagery.  These  are  generated  in 
cloudy  areas  and  are  detectable  in  the  NIR  due  to  the  shift  to 
smaller  size  cloud  droplet  distribution  from  the  addition  of 
CCN.  These  are  not  seen  in  the  visible  because  at  0.63  /^m 
there  is  no  absorption.   Therefore,  the  reflectance  becomes 


primarily  a  function  of  the  liquid  water  property  and  size 
distribution.  At  3.7  /um,  there  is  moderate  absorption  hence 
reflectance  becomes  primarily  dependent  on  the  particle 
radius.  At  11  /xm,  there  is  high  absorption  and  the  particle 
acts  as  a  black  body  with  no  reflectance. 

Fett  et  al.  (1979)  indicated  that  the  most  common  areas  of 
ship  track  formation  were  in  regions  of  closed  cellular  clouds 
with  medium  base  heights.  Morehead  (1988)  pointed  out  that 
other  conditions  which  help  to  promote  ship  track  formation 
include  areas  with  saturated  or  supersaturated  air  near  the 
top  of  the  marine  layer;  areas  where  the  air  temperature  is 
cooler  than  the  sea-surface  temperature;  or  areas  where  there 
is  a  minimum  of  vertical  wind  shear.  Saturated  air  at  the  top 
of  the  boundary  layer  provides  sufficient  moisture  to  form 
cloud  droplets  in  the  presence  of  the  increase  of  CCN  provided 
by  the  ship  exhaust.  The  warm  sea  surface  temperature 
enhances  track  formation  in  the  presence  of  a  moist  boundary 
layer  by  creating  an  environment  which  is  more  favorable  to 
fog  and  haze.  The  effect  of  wind  shear  in  the  boundary  layer 
would  be  to  disperse  the  CCN  provided  by  the  exhaust  resulting 
in  poorly  defined  or  nonexistent  cloud  lines.  Ship  track 
cloud  formation  will  be  further  enhanced  in  areas  of 
stratus/ stratocumulus  cloud  topped  boundary  layers,  eastern 
oceans,  and  high  latitudes. 


E.   OBJECTIVES  AND  ORGANIZATION 

The  goal  of  this  thesis  is  to  analyze  ship  track  data  over 
an  extensive  data  base  of  AVHRR  satellite  imagery  using  the 
ship  track  detection  algorithm  developed  by  Coakley  et  al. 
(1987)  and  modified  by  Morehead  (1988) .  This  study  will  have 
two  basic  objectives.  One  is  to  develop  a  summary  of  where 
geographically  the  ship  tracks  most  often  occurred  during  the 
study  period.  An  additional  focus  is  to  conduct  a  statistical 
evaluation  of  confirmed  ship  tracks  to  compare  the 
relationship  and  contrasts  between  AVHRR  channel  1 
reflectivity  and  AVHRR  channels  3  and  4  radiance  for  ship 
track  areas  versus  non-ship  track  areas. 

Chapter  II  will  discuss  the  data  used  in  the  study  as  well 
as  the  satellite  and  sensor  used  to  obtain  the  data.  Chapter 
II  will  discuss  the  analysis  used  in  the  two  parts  of  this 
study,  including  the  techniques  used  to  manipulate  the  data 
for  determining  ship  track  geographical  locations  and  the 
techniques  used  in  the  statistical  study.  The  results  of  the 
analysis  will  be  presented  in  Chapter  III  and  Chapter  IV  will 
contain  the  conclusions  which  are  drawn  from  the  results  and 
recommendations  for  future  work. 
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H.   DATA  PROCESSING  AND  ANALYSIS 

The  first  part  of  Chapter  II  describes  the  satellite 
platform  and  sensor  used  to  collect  the  data  for  this  study. 
Then,  the  source  and  format  of  the  data  will  be  discussed  as 
will  the  basic  signal  processing  and  calibration  which  is 
conducted  prior  to  use  by  the  ship  track  cloud  retrieval 
algorithm. 

A.  SATELLITE 

The  platforms  providing  the  data  for  this  study  were  the 
NOAA  9  and  10  polar  orbiting  satellite  series  sponsored  by  the 
National  Oceanic  and  Atmospheric  Administration  (NOAA) .  The 
satellites  fly  sun-synchronous  polar  orbits  at  an  altitude  of 
about  525  nautical  miles.  The  NOAA  9/10  series  provided 
coverage  for  different  times  of  the  day.  NOAA-9  satellite 
passes  were  generally  around  2300  UTC  (1500  local  standard 
time)  and  are  roughly  parallel  to  the  North  American  coast. 
NOAA-10  passes  were  at  1600  UTC  (0800  local  standard  time)  and 
cut  NE-SW  across  the  central  California  coast.  The  data  were 
collected  by  Scripps  Satellite  Oceanography  facility  in  La 
Jolla,  California. 

B.  SENSOR 

Radiance  data  on  the  NOAA  9/10  satellites  are  collected  by 
the   Advanced   Very   High   Resolution   Radiometer   (AVHRR) 


instrument.  This  sensor  provided  a  nadir  resolution  of  1km  by 
lkm  at  the  polar  orbit  altitude  of  525  N  mi. 

Table  1  shows  the  five  wave  bands  recorded  by  the  AVHRR 
sensor.   Channel  1  is  in  the  visible,  channels  2  and  3  are  in 


TABLE  1 

WAVEBANDS  OF  THE  ADVANCED 

VERY  HIGH  RESOLUTION  RADIOMETER  (AVHRR) 


Channel 

Wavebands  (/xm) 

Center 
Frequency  (/im) 

1 

0.58  -   0.68 

0.63 

2 

0.73  -   1.10 

0.83 

3 

3.55  -   3.93 

3.7 

4 

10.30  -  11.30 

11.0 

5 

11.50  -  12.50 

12.0 

the  near- infrared,  and  channels  4  and  5  are  in  the  thermal 
infrared  portion  of  the  electromagnetic  spectrum. 

Channel  2  was  used  to  prepare  the  satellite  pass  overview 
that  distinguished  what  areas  to  investigate  for  possible  ship 
tracks.  Channels  1,  3  and  4  were  used  in  the  ship  track 
detection  algorithm. 


C.   SIGNAL  PROCESSING  AND  CALIBRATION 

The  development  of  cloud  analysis  algorithms  as  part  of 
the  International  Satellite  Cloud  Climatology  Project  (ISCCP) 
motivated  the  collection  of  the  data  used  in  this  study.  All 
the  data  were  part  of  the  First  ISCCP  Regional  Experiment 
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(FIRE)  conducted  in  the  summer  of  1987.  All  the  high  density 
tapes  used  for  the  analysis  had  been  previously  navigated  for 
the  FIRE  study. 

The  original  data  came  on  high  density  tapes  in  the  form 
of  counts  of  radiant  energy  in  each  of  the  five  channels 
listed  above.  The  radiant  energy  data  were  calibrated  and 
converted  to  geophysical  parameters.  The  calibration  of 
channels  1,  3  and  4  is  discussed  briefly  below  and  is 
discussed  in  full  by  Allen  (1987) . 

1.  Channel  1 

This  is  the  visible  channel,  which  was  used  in  this 
study  to  search  for  areas  where  the  potential  for  finding  ship 
tracks  was  greatest.  Since  it  is  known  that  ship  tracks  are 
not  likely  to  be  found  in  cloud  free  regions  in  the  visible, 
these  areas  could  be  avoided.  The  calibration  of  channel  1 
data  is  accomplished  assuming  a  linear  relationship  between 
counts  received  by  the  sensor  and  the  reflectance.  The 
calibration  is  done  in  terms  of  albedo  and  the  results  are  in 
units  of  percent  of  reflectance. 

2.  Channel  4 

Channel  4  imagery  is  used  to  help  classify  cloud 
masses.  It  is  used  in  the  Coakley  algorithm  to  help  locate 
non-ship  track  pixels  with  similar  microphysical  properties  as 
nearby  ship  track  pixels.   Channel  4  data  result  mainly  from 
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thermal  emission.    The  data  are  converted  to  radiance 
measurements   with  units   of   W/m2-sr-cm_1   using   a   linear 
correlation  to  counts  of  radiant  energy. 
3 .   Channel  3 

The  channel  3  near-infrared  band  data  contains  the  most 
important  information  for  ship  track  detection.  Ship  tracks 
show  up  very  clearly  in  channel  3  while  they  may  not  be 
detected  in  other  wavelengths.  Channel  3  data  is  comprised  of 
both  reflectance  and  thermal  emittance  information  for 
daylight  passes.  As  discussed  by  Morehead  (1988),  two  types 
of  channel  3  data  are  used: 

a.  Channel   3  Radiance 

The  first  is  a  channel  3  radiance  which  contains 
both  the  reflectance  and  emission  contributions,  determined  by 
a  linear  correlation  between  counts  measured  and  radiance, 
with  units  of  W/m2-sr-cm"1. 

b.  Channel   3  Reflectance 

The  second  is  a  channel  3  reflectance  signature 
estimating  only  that  portion  of  the  measured  irradiance 
resulting  from  reflectance. 

The  reflectance  is  derived  from  the  channel  3 
radiance  and  the  channel  4  radiance  which  represents  the 
thermal  emission.  The  channel  4  data  are  used  to  estimate 
emission  in  channel  3.  This  emission  is  subtracted  from  the 
total  radiance  leaving  only  a  reflectance. 
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D.  SHIP  TRACK  GEOGRAPHICAL  STUDY 

One  focus  of  this  analysis  was  to  determine  geographically 
where  the  ship  tracks  formed  during  the  period  of  the  study 
(1-19  July,  1987)  .  Each  AVHRR  data  tape  processed  contained 
an  overview  image  consisting  of  one  satellite  pass  of  the 
North  Pacific  Ocean  basin  (Figure  1)  .  Each  pass  was 
subjectively  scanned,  subscene  by  subscene  (500km  by  500km) , 
in  both  channel  3  and  channel  1  to  locate  ship  track  clouds 
(Figures  2  and  3) .  Once  a  ship  track  was  found,  latitude  and 
longitude  positions  were  traced  along  the  ship  track.  The 
ship  tracks  were  then  plotted  to  show  their  location. 
Geographical  diagrams  were  produced  for  each  pass  showing  the 
position  and  track  of  all  ship  tracks  found  in  both  channel  1 
and  channel  3 .  Figure  4  and  Figure  5  show  the  geographical 
plots  of  ship  tracks  found  in  channel  1  and  channel  3 
respectively,  for  the  7  July  pass  whose  overview  is  shown  in 
Figure  1. 

E.  STATISTICAL  STUDY 

The  second  part  of  this  study  involved  a  statistical 
comparison  of  channels  1,  3  and  4.  An  algorithm  developed  by 
Coakley  et  al.  (1987)  and  modified  by  Morehead  (1988)  was 
applied  to  each  subscene  determined  to  contain  ship  tracks, as 
found  in  the  geographical  process  above.  The  result  of  each 
subscene  statistical  analysis  was  compared  to  the  respective 
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Figure  1.  Channel  2  overview: 
NOAA-9  satellite  pass  AR6085  (2237 
UTC,  7  July  1987)  .  Displays  the 
location  of  the  subscene  used  as  an 
example  of  the  data  processing  of 
this  thesis. 
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Figure  2.  Channel  3  near-infrared  subscene:  AR6085  (2237  UTC, 
7  July  1987) . 
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Figure  3.   Channel  1  visible  subscene: 

7  July  1987) . 


AR6085  (2237  UTC, 
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Figure  4.   Channel  1  ship  tracks:   AR6085  (22  37  UTC, 
7  July  1987) . 
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Figure  5.   Channel  3  ship  tracks:   AR6085  (22  37  UTC, 
7  July  1987) . 
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channel  3  image.  A  numbered  8  by  8  grid  of  64  pixel  squares 
was  overlaid  on  both  the  statistical  image  and  the  channel  3 
image.  The  numbered  squares  (1  -  64)  facilitated  in  selecting 
which  pixel  areas  to  filter  out.  Only  the  grid  areas  with 
confirmed  ship  track  data  were  reprocessed,  filtering  out  non- 
ship  track  statistical  data. 

This  process  is  shown  for  the  NOAA-9  pass  on  7  July  1987, 
in  Figures  6,  7  and  8.  Figure  6  shows  a  channel  3  subscene 
with  ship  tracks  clearly  visible  as  dark  curvilinear  lines 
with  an  8  by  8  grid  overlaid.  Figure  7  shows  the  results  of 
the  Coakley/Morehead  algorithm  processed  on  the  subscene. 
Figure  8  shows  the  reprocessed  analysis  with  the  non-ship 
track  lines  filtered  out  and  the  normal  or  ambient  control 
points  included.  The  filtering  process  was  conducted  to 
create  a  statistical  data  base  representing  as  near  as 
possible,  only  confirmed  ship  tracks.  The  final  step  in  the 
analysis  was  to  process  grid  points  containing  confirmed  ship 
tracks  through  a  comparison  routine  of  ship  track  vs.  ambient 
pixels,  giving  values  of  reflectance  for  channel  1  and 
radiance  for  channels  3  and  4.  The  results  of  this  routine 
were  presented  in  tabular  and  scatter  plot  format  with 
comparisons  among  channels  1,  3  and  4.  Figures  9  and  10  are 
the  scatter  plots  for  the  case  presented  in  Figures  6  through 
8  above.  Figure  9  shows  the  comparison  of  channel  1 
reflectivity  to  channel  3  radiance.  Ambient  pixels  are 
represented  by   (*)s  an<*  enclosed  by  a  solid  ellipse. 
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Ship  track  pixels  are  represented  by  (+)s  and  enclosed  by  a 
dashed  ellipse.  The  increase  in  reflectivity  in  channel  1  and 
radiance  for  channel  3  for  the  ship  track  pixels  can  be 
readily  seen  in  Figure  9.  Figure  10  shows  a  comparison 
between  channel  1  reflectivity  and  channel  4  radiance.  It  is 
noted  that  there  is  very  little  difference  between  the  ambient 
and  ship  track  pixel  values  in  this  comparison.  This  result 
will  be  discussed  further  in  Chapter  III. 
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Figure  6.  Channel  3  subscene  with  8X8  grid  overlay:  AR6085 
(2237  UTC,  7  July  1987) . 


21 


1  M 

.. . 

T 

rj  ? 

P 

a 

64 

63 

62 

V  61 

60 

^ 

58 

57 

'■;■■■  ? 

56 

V 

55  *- 

t 

"SS 

4" 

■r    52 

■■  ,y  51 . 

^0 

49 

> 

48 

47 

445 

44 

i'  '■, ': 

^  43 

42  i 

41 

,1 

IT 

f 

i   40 

39 

38 

37 

t   36; 

35 

J     34 

33 

7 

* 

31 

~4 

r 

32 

30 

29 

28^ 

►'  27 

26 

25 

flP 

a 

24 

23 

22 

21  < 

20 

IS 

IB  ; 

17 

16 

15 

14 

13 

12 

11 

10 

9 

8 

7 

6 

5 

4 

3 

2 

1 

Figure  7.   Coakley/Morehead  algorithm  on  Fig  6.  subscene. 
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Figure  8.   Fig.  7  subscene  reprocessed:   After  grids 

containing  non-ship  track  analysis  were 
filtered  out. 
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Figure  9.   Scatter  plot,  chn  1  reflectivity  vs.  Chn  3 

radiance:   Ship  track  pixels  (+)  and  ambient 
pixels  ( • ) . 
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Figure  10.   Scatter  plot  Ctan  1  reflectivity  vs.  Chn  4 

radiance:   Ship  track  pixels  (+)  and  ambient 
pixels  ( • ) . 
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m.  RESULTS 

A.   SHIP  TRACK  FREQUENCY  AND  GEOGRAPHICAL  ANALYSIS 

Coakley  et  al.  (1987)  made  an  assessment  that  ship  tracks 
would  occur  in  only  5  to  10%  of  all  orbital  passes.  During 
the  course  of  this  study,  28  NOAA  passes  (17  NOAA  -  9, 
11  NOAA  -  10)  were  analyzed.  This  encompassed  17  days  from  1 
to  19  July  1987.  During  this  period,  408  subscenes  (500km  by 
500km)  were  examined  leading  to  the  discovery  of  316  ship 
tracks.    Table  2  provides  a  summary  of  the  ship  tracks  found 


TABLE  2 
RESULTS  OF  GEOGRAPHICAL  ANALYSIS 

PASSES  WITH  SHIP  TRACKS 

SATELLITE 

PASSES 

IN  CHN  3 
ONLY 

IN  CHNs 
1  &  3 

TOTAL 

NOAA- 9 

17 

3 

10 

13  (76%) 

NOAA- 10 

11 

5 

4 

9  (82%) 

TOTAL 

28 

8  (29%) 

14  (50%) 

22  (79%) 

in  the  analysis.  Only  six  of  the  passes  did  not  contain  ship 
tracks  in  either  channel  1  or  3.  This  analysis  found  that  79% 
of  all  passes  studied  contained  ship  tracks,   which  is 
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Figure  11.   Channel  1  ship  tracks,  1-19  July  1987. 

dramatically  higher  than  the  Coakley  et  al.  assessment. 
Figures  11,  12  and  13  display  the  plots  of  ship  tracks  found 
in  channel  1,  3  and  combined  respectively  for  all  the  passes 
processed  for  this  study. 
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Figure  12.   Channel  3  ship  tracks,  1-19  July  1987 


28 


Figure  13.   Chnl/Chn3  ship  tracks  combined,  1-19  July 
1987. 
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B.   STATISTICAL  STUDY 

Ship  track  studies  by  Coakley  et  al.  (1987)  and  Morehead 
(1988)  were  conducted  on  special  cases  with  only  a  limited 
data  set  and  provided  a  snapshot  of  the  radiative  properties 
of  a  few  ship  tracks.  The  analysis  of  this  thesis  was 
conducted  over  a  much  broader  data  base.  As  the  geographical 
analysis  was  being  conducted,  all  subscenes  found  to  contain 
ship  tracks  were  annotated  for  further  analysis  using  the 
modified  Coakley  algorithm.  Once  the  filtering  process, 
discussed  previously,  was  completed,  a  statistical  data  base 
was  created  that  provided  a  data  enriched  view  of  the 
radiative  properties  of  the  ship  tracks  detected. 

One  of  the  results  from  the  algorithm  was  the  degree  of 
variability  in  the  radiative  and  reflective  properties  of  the 
ship  tracks.  Table  3  presents  the  maximum  and  minimum  values 
of  channel  1  reflectivity  and  channels  3  and  4  radiance  for 


TABLE  3 

AMBIENT 
MAX 

PIXELS 
MIN 

SHIP  TRACK  PIXELS 
MAX        MIN 

CHANNEL  1 

39.5942 

16.9879 

45.1522 

21.2278 

CHANNEL  3 

0.9252 

0.4462 

1.0548 

0.5400 

CHANNEL  4 

96.6457 

76.6413 

96.6786 

76.6425 
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the  ambient  and  ship  track  pixels  processed.  Table  4  gives 
mean  radiance  values  for  the  channel  1,  3  and  4  ambient  and 
ship  track  pixels,  for  the  61  subscenes  processed.  When  the 
maximum  and  minimum  values  shown  in  Table  3  are  compared  to 
the  average  of  the  means  from  Table  4 ,  one  can  see  that  for 
the  ship  track  pixels,  the  spread  from  the  average  is  greater 
than  the  spread  from  the  average  for  the  ambient  pixels.  This 
implies  that  the  response  of  clouds  to  ship  effects  is  not 
confined  to  one  resulting  reflectance  value. 

Figures  14  through  35  represent  some  selected  scatter 
plots  from  the  data  set.  The  diagrams  plot  distributions  of 
channel  1  reflectivity  as  a  function  of  channel  3  or  channel 
4  radiance.  Pixels  contaminated  by  ship-stack  exhaust  are 
designated  by  (+)  and  those  of  randomly  selected,  nearby 
ambient  fields  of  view,  are  designated  by  (•)•  Tne  solid 
(ambient  control  group)  and  dashed  (ship  track)  ellipses 
represent  fits  at  2  standard  deviations  to  the  reflectivity 
and  radiance  distributions,  which  are  taken  to  be  joint 
Gaussian  distributions. 

The  study  by  Coakley  et  al.  (1987)  displayed  a  similar 
scatter  plot  depicting  an  overall  increase  in  both  channel  1 
reflectivity  and  channel  3  radiance  for  the  selected  subscene 
of  that  study.  It  can  be  seen  in  Figures  14  through  35  that 
the  visible  reflectance  is  not  always  increased  for  the  ship 
track  group.  In  fact,  only  36  of  the  61  files  presented  in 
Table  4,  had  increases  in  channel  1  reflectivity.  Figures  15, 
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TABLE    4 

HI.E  AMHIKIIT    I'lXKIi    MEANS  SHIP    TRACK    PIXEL    MEANS 

Cllll    i  cim   3  ciih  4  cum  aim  cim4 
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.1293 

26 
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0.7312 

88 

.0581 

9 
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.  5255 

0. 6044 

85 

.7383 

26 

.4225 

0. 7  2  69 

85 

.8990 

10 

21 

.5716 

0.6344 

81 

.7598 

24 

.8219 

0.8223 

81 

.9286 

11 

23 

.  6894 

0. 4869 

87 

.0027 

26 

.5680 

0.6166 

86 

.9503 

12 

19 
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0. 5098 

87 

.2836 

23 

.5276 

0.7478 

87 
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21 
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94 

.4822 

21 

.9398 

0.6317 

94 

.5021 
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21  -  25,  and  27  all  show  a  decrease  for  the  visible 
reflectance.  The  scatter  plots  graphically  show  the  great 
deal  of  variability  in  both  the  ship  track  and  control  group 
values.  For  example,  Figure  16  shows  a  range  in  channel  1 
ambient  reflectance  of  12  -  32%  and  in  channel  3  ambient 
radiance  of  0.5  -  0.75  mWm"2sr'1cm.  Figure  16  shows  a  range  in 
channel  1  ship  track  reflectance  of  19  -  31%  and  in  channel  3 
ship  track  reflectance  of  0.71  -  0.94  mWm^sr'cm.  Coakley 
et  al.,  attributed  this  variability  to  the  irregularity  in  the 
distribution  and  size  of  cloud  droplets  within  the  same  cloud 
mass,  which  serves  to  obscure  some  of  the  increase  in  channel 
1  reflectivity.  It  is  also  shown  in  Figures  14  through  35, 
that  the  channel  3  near-infrared  radiance  values  are  higher  in 
every  case  represented  in  this  study,  confirming  similar 
observations  made  by  Coakley  et  al. 

Table  4  shows  the  small  variability  of  the  channel  4 
radiance  and  also  the  almost  negligible  variation  between  ship 
track  and  control  group  radiance  levels.  Coakley  et  al., 
attributed  this  to  the  observation  that  at  the  channel  4 
wavelength,  liquid  water  is  a  strong  absorber,  therefore  the 
amount  of  scattering  is  negligibly  small.  Also,  there  are 
enough  droplets  so  that  marine  stratocumulus  emit  like  black 
bodies  regardless  of  droplet  size  and  concentrations.  Thus, 
there  is  no  change  in  thermal  emission  between  contaminated 
and  noncontaminated  clouds. 
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Figure   14.      Scatter  plot  from  Table  4   -  File 

number    8.       AR6105    (2314    UTC,    13/7/87). 
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Figure  15.   Scatter  plot  from  Table  4  -  File 

number  9.  AR6105  (2314  UTC,  13/7/87). 
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Figure  16.   Scatter  plot  from  Table  4  -  File 

number  10.   AR6105  (2314  UTC,  13/7/87) 
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Figure  17.   Scatter  plot  from  Table  4  -  File 

number  11.   AR6105  (2314  UTC,  13/7/87). 
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Figure   18.      Scatter  plot   from  Table  4   -  File 

number  12.   AR6105  (2314  UTC,  13/7/87) 
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Figure  19.   Scatter  plot  from  Table  4  -  File 

number  13.   AR6105  (2314  UTC,  13/7/87). 
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Figure  20.   Scatter  plot  from  Table  4  -  File 

number  14.   AR6113  (1526  UTC,  15/7/87). 
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Figure  21.   Scatter  plot  from  Table  4  -  File 

number  15.   AR6114  (2252  UTC,  15/7/87). 
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Figure  22.   Scatter  plot  from  Table  4  -  File 

number  16.   AR6114  (2252  UTC,  15/7/87). 
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Figure  23.   Scatter  plot  from  Table  4  -  File 

number  17.   AR6114  (2252  UTC,  15/7/87). 
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Figure  24.   Scatter  plot  from  Table  4  -  File 

number  18.   AR6120  (2231  UTC,  17/7/87). 
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Figure  25.   Scatter  plot  from  Table  4  -  File 

number  19.   AR6120  (2231  UTC,  17/7/87). 
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Figure  26.   Scatter  plot  from  Table  4  -  File 

number  20.   AR6049  (1529  UTC,  1/7/87). 
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Figure  27.   Scatter  plot  from  Table  4  -  File 

number  21.   AR6049  (1529  UTC,  1/7/87) . 
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Figure  28.   Scatter  plot  from  Table  4  -  File 

number  22.   AR6049  (1529  UTC,  1/7/87). 
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Figure  29.   Scatter  plot  from  Table  4  -  File 

number  23.   AR6049  (1529  UTC,  1/7/87). 
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Figure  30.   Scatter  plot  from  Table  4  -  File 

number  24.   AR6049  (1529  UTC,  1/7/87). 
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number  25.   AR6051  (2213  UTC,  1/7/87). 
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Figure  32.   Scatter  plot  from  Table  4  -  File 

number  26.   AR6051  (2213  UTC,  1/7/87) . 
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Figure  33.   Scatter  plot  from  Table  4  -  File 

number  27.   AR6054  (1508  UTC,  2/7/87). 
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Figure  34.   Scatter  plot  from  Table  4  -  File 

number  28.   AR6054  (1508  UTC,  2/7/87). 
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Table  5  displays  the  differences  in  the  ship  track  pixel 
means  from  the  ambient  pixel  means.  Table  5  also  shows  the 
conversion  of  channel  3  radiance  to  reflectivity.  Table  5 
shows  that  there  is  a  marked  increase  in  values  for  channels 
1  and  3  in  the  ship  track  pixels  compared  to  ambient  pixels. 
Channel  1  ship  track  differences  increased  by  1.3%  compared  to 
ambient  pixel  means.  Channel  3  ship  track  differences 
increased  by  18.6%  compared  to  the  ambient  pixels.  One  can 
also  see  that  there  is  only  a  very  slight  increase  in 
channel  4.  This  is  again  in  support  of  the  findings  of 
Coakley  et  al.  The  percentage  of  increase  in  the  values  of 
reflectivity  and  radiance  was  found  to  be  less  in  this  study 
than  that  found  in  the  results  from  Coakley  et  al.,  as  seen  in 
Table  6.  Coakley  et  al.,  assumed  that  comparing  the  ratio  of 
change  in  reflectivity  for  channel  1  to  that  for  channel  3, 
removed  the  effect  of  fractional  cloud  cover  which  could 
influence  the  changes  in  the  reflectivities  noted  above.  The 
ratio  of  change  in  reflectivity  for  this  analysis  was 
0.18  +/-1.4  which  is  nearly  half  the  value  of  0.4  +/-  0.8 
found  by  Coakley  et  al. 
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TABLE  6 
COMPARISON  TO  COAKLEY  et  al.  (1987) 

CHANGE  IN 

RADIANCE 

(mWm^-sr'-cm) 

% 

REFLECTIVITY 

INCREASE 

RATIO  OF 
CHANGE  IN 
REFLECTIVITY 

CHANNEL  4 

CHANNEL  1 

CHANNEL  3 

CHN1  /  CHN3 

COAKLEY  et  al.  (1987) 

0.0+/-0.05 

1.6+/-0.7 

3.9+/-0.4 

0.4+/-0.8 

LUTZ 

0.0+/-0.12 

0.34+/-2.6 

2.1+/-0.99 

0.16+/-2.6 
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IV.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  primary  focus  of  this  thesis  involved  evaluating  a 
large  satellite  data  set  to  determine  the  number  and 
geographical  position  of  ship  tracks  that  occurred  during  the 
period  of  study  (1-19  July  1987) .  In  total,  316  ship  tracks 
were  located  and  plotted.  A  statistical  evaluation  was  then 
conducted  on  the  radiative  properties  of  the  confirmed  ship 
tracks . 

It  was  discovered  that  22  of  the  28  (76%)  satellite  passes 
evaluated  in  this  thesis  contained  ship  tracks.  The  76%  was 
shown  to  be  significantly  higher  than  the  5  to  10%  expected  in 
the  earlier  by  Coakley  et  al.  (1987) .  It  was  shown  that  the 
AVHRR  channel  3  was  more  favorable  to  the  observance  of  ship 
tracks  than  channel  1. 

The  large  number  of  ship  tracks  found,  provided  a  rich 
data  set  which  was  evaluated  using  the  modified  Coakley 
algorithm  for  automatic  ship  track  detection.  The  statistical 
study  concurred  with  the  conclusions  of  Coakley  et  al.,  that 
ship  exhaust  increases  the  reflective  properties  of  clouds. 
Table  4  indicated  there  is  substantial  variability  of  the 
cloud  radiative  properties,  which  show  up  clearly  when 
examining  the  ship  track  and  ambient  pixels  for 
channels  1,  3  and  4.  When  considering  the  means,  all  three 
channels  showed  an  increase  in  the  ship  pixels  compared  to  the 
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ambient  pixels.  The  increase  for  channels  1  and  3  was 
significant  (0.35  +/-  2.6%  reflectivity  for  channel  1  and 
0.11  +/-  0.05  mWm^sr'cm  for  channel  3)  and  very  slight  for 
channel  4  (0.01  +/-  0.12  mWm^sr^cm)  .  Although  all  three 
channels  showed  increases  in  radiative  properties,  as  in 
Coakley  et  al.,  the  increases  were  less  for  this  study. 

Many  areas  of  research  have  yet  to  be  explored  in  the 
study  of  ship  tracks.   Some  of  the  areas  for  future  work  are: 

1.  Investigation  of  weather  patterns  within  areas  of  known 
ship  track  formation  should  be  conducted  to  determine 
more  precisely  what  factors  encourage  the  formation  of 
ship  tracks. 

2 .  An  effort  should  be  made  to  discover  the  cause  of  the 
variability  found  in  the  reflectivity  and  radiance 
readings  associated  with  the  ship  track. 

3.  Work  should  be  continued  to  improve  the  ability  of  the 
Coakley  algorithm  to  process  satellite  data  and  auto- 
matically locate  and  analyze  ship  tracks.  Additionally 
the  algorithm  should  be  developed  to  automatically  create 
a  plot  of  ship  tracks  found. 

4.  Studies  are  still  needed  to  determine  if  specific  ship 
types  can  be  associated  with  variations  in  ship  track 
properties. 

The  study  of  ship  track  formation  will  be  a  continuing 
source  of  exciting  research  in  the  years  to  come.  Civilian 
and  military  use  of  the  ability  to  monitor  the  movements  of 
ships  by  evaluating  their  ship  tracks  is  very  promising.  With 
the  growing  worldwide  concern  about  our  environment,  the 
efforts  expended  on  ship  track  formation  processes  will  pay 
great  dividends  in  our  efforts  to  protect  our  fragile  planet. 
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